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In purely repulsive, C41, -symmetric Hubbard clusters a correlation effect produces an effective two- 
body attraction and pairing; the key ingredient is the availability oiW = Q pairs, that is, two-body 
solutions of appropriate symmetry. We study the tunneling of bound pairs in rings of 5-site units 
connected by weak intercell links; each unit has the topology of a Cu04 cluster and a repulsive 
interaction is included on every site. Further, we test the superconducting nature of the response 
of this model to a threading magnetic field. We present a detailed numerical study of the two-unit 
ring filled with 6 particles and the three-unit ring with 8 particles; in both cases a lower filling yields 
normal behavior. In previous studies on Id Hubbard chains, level crossings were reported (half- 
integer or fractional Aharonov-Bohm effect) which however cannot be due to superconducting pairs. 
In contrast, the nontrivial basis of clusters carrying VK = pairs leads to genuine Superconducting 
Flux Quantization (SFQ). The data are understood in terms of a cell-perturbation theory scheme 
which is very accurate for weak links. This low-energy approach leads to an effective hard core 
boson Hamiltonian which naturally describes itinerant pairs and SFQ in mesoscopic rings. For the 
numerical calculations, we take advantage of a recently proposed exact diagonalization technique 
which can be generally applied to many-fermion problems and drastically reduces the size of the 
matrices to be handled. 

71. 27. -fa Strongly correlated electron systems; heavy fermions 
74.20.Mn Nonconventional mechanisms 

73.22.-f Electronic structure of nanoscale materials: clusters, nanoparticles, nanotubes, and 
nanocrystals 



I. INTRODUCTION 

Low-dimensional systems with strong electron-electron cor- 
relations may lead to an anomalous Aharonov-Bohm (AB) 
effect [1] with ground-state energy oscillations versus flux <j) 
having a period shorter than the foundamental one, given by 
4>o = hc/e. However, as we shall see below, a fractional pe- 
riod (jto/N does not generally mean that the current is carried 
by particles with an effective charge e* — Ne. In particular, 
if A'^ = 2 the half-integer AB effect is not equivalent to the 
superconducting flux quantization (SFQ) and pairing is not 
necessarily implied. 

The repulsive Hubbard ring in the presence of a magnetic 
flux has been studied by several authors. At half filling it 
was inferred by numerical evidence [2] that the period of the 
ground state energy (0) as a function of the trapped flux 
(j> is a whole fluxon if U /t is in a physical parameter range. 
Away from half filling Kusmartsev [3] pointed out that for 
microscopic Hubbard rings with N particles the ground state 
energy has N minima in the range [0, ^o) for U = 00 and at 
low density. [This result has been confirmed by Schofield et 
al. in Ref. [4] .] These oscillations were originally explained in 
terms of spin-fiip processes by the same author, so that the 
system varies its total spin as the fiux changes by one fiux 
quantum. However, when U — 00 many degenerate ground 
states with different total spin S exist, and a more accurate 
explanation of the fractional AB effect has been provided by 
Yu and Fowler [5] [6] . They studied the Lieb- Wu [7] equations 
for a chain with twisted boundary conditions [a fiux (j> corre- 
sponds to a twist of 2TT<l)/(j}o] with leading t/U corrections. 
Increasing the flux the holon momenta get shifted and the 



energy of the holon sea grows; to counterbalance this effect 
the system generates a compensating momentum by creat- 
ing a hole in the distribution of the spinon quantum num- 
bers. This excitation in the spinon sea is not energetically 
suppressed as far as Nt/LU <^ 1 where L is the number of 
sites. The quantization of the spinon momenta does not al- 
low a full compensation for the effect of a continuosly varying 
flux (j> and leads to the energy oscillations discussed above. 
As Nt/LU increases the spinon excitation energies are raised 
and beyond a critical value Eo{(t>) has only two minima at 
(j) — and = (^0/2 corresponding to the spinon excitations 
with momentum and tt, as numerically observed in Refs. [2] 
[5] [8]. Hence, the half-integer AB effect is driven by level- 
crossing due to the spinon degrees of freedom and it is not 
SFQ. 

Recently, an exact result on the inflnite U Hubbard ring 
with twisted boundary conditions has been obtained by 
Nakano [9]. He proved that the ground state energy E^'^\(j>) 
with an even number of particles N is periodic with period 
(j)o/2 if the z-component of the total spin Sz = and with 
period (po/N if M|/M| (with M| > Af|) is not an integer; 
here Ma is the number of particles with spin a =t,i- We 
emphasize that the Nakano theorem is not in contrast with 
the results by Kusmartsev. Indeed, even ff £;(")((?!)) has period 
0o/2 it may exhibit other minima in the range [0, 0o). 

The half-integer AB effect caused by the existence of con- 
flned pairs has been observed in the framework of extended Id 
Hubbard models. Sudb0 et al. [10] studied an alternating Cu- 
O ring with a charge-gap between Cu and O (which favours 
the particles to occupy O sites) and an on-site repulsion is 
allowed on Cu's but not on O's. In this way a half-integer 
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flux quantization needs a strong ofT-site repulsion; pair con- 
finement is achieved by allowing for repulsion-free sites, and 
not by producing an effective attraction. Other authors con- 
sidered modified Id models with bond-charge interactions so 
that the number of doubly occupied sites is an extra con- 
served quantity [11]. Turning on a term breaking the extra 
symmetry gives rise to an effective attractive interaction and 
to the SFQ for moderate on-site repulsion. However, the pair- 
ing mechanism is not driven by the Hubbard interaction as 
witnessed by the presence of SFQ even at ?7 = [12] [13]. 

Few data are available in 2d Hubbard systems, since no 
exact solutions exist in this case and numerical explorations 
are possible only for small clusters. Long ago Canright and 
Girvin [14] discussed the magnetic response of the Hubbard 
model showing SFQ in the attractive case by threading a 
magnetic field into a cylindrical probe. Assaad et al. [15] 
computed numerically the magnetic response of the repulsive 
Hubbard model in the same geometry at quarter filling, but 
the flux quantization was found to be normal. Signatures of 
anomalous flux quantization in a 4 x 4 geometry have been 
provided by Arrachea et al. [16] ; however they studied an ex- 
tended Hubbard model with nearest-neighbour correlated, e.g. 
occupation-dependent, hopping. 

All these results would suggest that the half-integer AB ef- 
fect cannot be interpreted as SFQ in the context of repulsive 
Hubbard systems. However, Anderson [17] first advocated 
the non-conventional superconductivity arising from repulsive 
interactions, proposing the one-band Hubbard Hamiltonian 
as a prototype model. Evidence for pairing has then been 
obtained by several authors by a variety of methods. Ana- 
lytic approaches based on a renormalization method [18] [19] 
and on various implementations of the renormalization group 
technique [20] [21] [22], generalized conserving approximation 
theories like FLEX [23], as well as Quantum Monte Carlo 
studies on supercells [24] lead to this conclusion. In Refs. 
[25] [26] [27] [28] we show that pairing occurs also in purely 
repulsive, C4i;-symmetric Hubbard clusters. 

In this paper we wish to study the tunneling of bound pairs 
in rings of 5-site units with a Cu04 topology, see Fig.(l). 




FIG. 1. Illustration of the ring topology described in the main 
text. T and tcu represent O-O and Cu-Cu links, respectively, as 
explained in detail in Section III below. 

In the following we shall refer to the central site as Cu and 



to the four external sites as O just to distinguish their posi- 
tion in the unit cell. We shall see that numerical solutions 
of such a model clearly show superconducting pair hopping 
if the total number of particles is 2|A| -I- 2p where |A| is the 
immber of units and < p < |A|; in particular, once a mag- 
netic field is switched on into the ring, SFQ is unambiguously 
observed. We note that in our model no superconducting re- 
sponse is obtained with less than 2 particles per 5-site unit 
[CUO4]. This means that the Zhang- Rice picture [29] for the 
two-dimensional d — p model does not work in the present 
context and we wish to explore a scenario with a larger num- 
ber of particles. The data are interpreted by implementing 
a cell-perturbation theory which is very accurate for weak 
links. This low-energy approach leads to an effective hard 
core boson Hamiltonian which naturally describes itinerant 
pairs and SFQ in mesoscopic rings as well. We feel that SFQ 
from purely repulsive Hubbard models is interesting by itself 
even if it may have no direct relevance for high- Tc cuprates. 
The tunneling of pairs has been studied in the context of 
the t-J model too [30], but here we wish to study symmetric 
Hubbard clusters for a broad range of U /t including the weak 
coupling limit. 

The plan of the paper is the following. In Section II we 
recall that the CUO4 cluster has two-body singlet eigenstates 
without double occupation called W — pairs. By exploiting 
a symmetry driven configuration interaction mechanism they 
get bound once dressed by the virtual electron-hole excita- 
tions; the binding is due to an effective attraction among the 
particles of the W = pair. In Section HI we introduce the 
model Hamiltonian with symmetric clusters, as units of an 
arbitrary graph, linked by intercell hoppings. For any given 
graph, we deduce a low-energy effective Hamiltonian in Sec- 
tion IV. Partial results presented elsewhere [31] are extended 
to any number of conducting pairs and we show that the ef- 
fective Hamiltonian is equivalent to the Id antiferromagnetic 
Heisenberg-Ising model for a ring topology. In Section V we 
recall a recently proposed spin-disentangled diagonalization 
technique that allowed us to carry out the numerical calcu- 
lations. We also give practical details on the performance of 
this novel approach. Results are presented and interpreted 
physically in Section VI for the two-unit ring and in Section 
VII for the three-unit ring and compared with the analytic 
predictions. For the larger ring, after reviewing from Ref. 
[31] the flux dependence of the energy levels in the case of 
O-O intercell bonds we analyze the current and show that 
it is a discontinuous function of the flux and can be param- 
agnetic as well as diamagnetic. Novel results for the case of 
Cu-Cu bonds are presented, and this case is definitely more 
exotic: the bound pairs have an infinite effective mass and 
the flux quantization is normal. Our main conclusions are 
summarized in Section VIII. 



II. W = PAIRING IN HUBBARD MODELS 

The repulsive Hubbard Hamiltonian has two-body singlet 
eigenstates without double occupation [25] [26] [27] [28] called 
W = pairs. Such solutions are also allowed in the fully 
symmetric clusters C. In the many-body ground state these 
pairs get dressed and bound, and this is signaled by Ac{N) < 
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where Ac{N) = E^^'>{N) + E^^'>{N - 2) - 2E^^'>{N - 1); 
E^-' (N) is the interacting ground state energy of the cluster 
C with N particles. By means of a non-perturbative canonical 
transformation [18] [19], it can also be shown that Ac{N) < 
is due to an attractive effective interaction and at weak 
coupling |Ac(A'^)| is just the binding energy of the pair. The 
extension of the theory to the full plane was also put forth in 
Ref. [18]. 

The C4v symmetric 5-site cluster is the smallest one where 
the W = pairing mechanism works. We have already de- 
scribed W — pairing in great detail as a function of the 
one-body and interaction parameters on all sites; the study 
was extended to larger clusters too [27] [32]. The 5-site unit 
has the same topology of the Cu04 cluster; thus, we label the 
central site by the Cu symbol and the four external ones by 
the O symbol. The 5-site cluster will be also called CUO4. 
In order to simplify the analytical formulas, we neglect the 
O-O hopping term and the only nonvanishing hopping ma- 
trix elements are those between an O site and the central Cu 
site; they are all equal to t. For the sake of simplicity, we 
parametrize the Hubbard model in such a way that actually 
everything depends only on the ratio U/t. Thus, we consider 
the Hubbard Hamiltonian 



HcuC 
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where pj^ and pi^ are the creation and annihilation opera- 
tors onto the O i — 1,..,4 with spin a =Ti J.) <ia do- 
are the creation and annihilation operators onto the Cu site, 
while h[^^ = Pia-P''^ ^^'^ '^o-'*' ~ d-lrdc are the corresponding 
number operators. -H'cu04 is invariant under the permuta- 
tion group 54. The classes of 54 can be labelled with the 
quatern (711,712,713,714). Each class contains all the permuta- 
tions where 71; is the maximum number of sets of i elements 
that remain unchanged after the permutation. For example 
(4,0,0,0) contains all the permutations such that 4 sets con- 
taining one element remain unchanged, that is the identity. 
On the other hand (2,1,0,0) contains all the permutations such 
that 2 sets containing one element and one set containing 2 
elements remain unchanged. Starting from the configuration 
1234, this class contains the permutations 1243, 1432, 1324, 
4231, 3214, 2134. S4 has the irreducible representations {ir- 
reps) Ai (total-symmetric), B2 (total-antisymmetric), £ (self- 
dual), Ti and its dual T2, of dimensions 1, 1, 2, 3 and 3, 
respectively, see Table I. The model admits a 14^ = pair 
belonging to the irrep £ and formed by mixing degenerate 
one-body states. 
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Table I: Character Table of the 54 group. 
(711,712,713,714) labels a class of S'4. 
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For later use, we recall how the irrep F G 
{Ai, T\, £, T2, B2} of S4 breaks in Civ, that is the point 
symmetry group of the square. Cav is a subgroup of S'4 and 
its Character Table is shown in Table II. From Table I and II 
we have 
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T2=A2®E, 62= B2, 
> B2 . (2) 
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Table II: Character table of the C4v symmetry group. Here 
1 denotes the identity, C2 the 180 degrees rotation, 0^"*"^ c| ' 
the counterclockwise and clockwise 90 degrees rotation, cfx , Oy the 
reflection with respect to the 7; = and a; = axis and a^, a— the 
reflection with respect to the x = y and x = —y diagonals. In the 
last column we show typical basis functions. 



/.From Eq.(2) we see that we may label the two components 
of the W — 0-pair irrep £ in terms of the irreps Ai and B2 
of C4v For example, defining by p^^^ = "Tf o- ~ P3 o-) ^^'^ 
p]^ o — -^(P2<T ~ p\ cr) ^'^o of the three degenerate eigen- 
operators of the kinetic term in HcuOi, the B2 component 
of the = pair is obtained by acting with the two-body 
singlet operator 
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on the vacuum. 

The ground state of Cu04[2] (i.e. CUO4 with 2 particles) 
belongs to ^^1 and that of Cu04[4] is in ^£; both are singlets, 
as the notation implies. The ground state of CUO4 [3] is a ^7i 
doublet. Acu04(4) has a minimum at i7 « 5 t for this model, 
as shown in Fig. (2), and it is negative when < U < 34.77 t. 
We emphasize that Acu04(4) becomes positive for large val- 
ues of U /t and hence pairing disappears in the strong coupling 
regime. In the present problem U must exceed several tens 
of times t before the asymptotic strong coupling regime sets 
in. A perturbation theory will strictly apply at weak coupling 
where the second derivative of the curve is negative; however, 
qualitatively a weak coupling approach is rewarding in all the 
physically interesting range of parameters. The above pairing 
mechanism does not work in the neighborhood of the infinite 
U limit. 
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FIG. 2. Acu04(4) (in t units) as a function of U/t. The 
maximum binding occurs at f ~ 5t where Ac^q^{4) « —0.042 t. 
For U > 34.77 t (not shown), Acu04(4) becomes positive and 
pairing disappears. 



the Hubbard model belongs to the broad cathcgory of Kohn- 
Luttinger effects with no direct interaction and a second-order 
correction coming from the spin-flip channel only, as shown 
in Fig. (3). 



We analyzed the pairing mechanism in detail in Ref. [28]; 
for the sake of simplicity hero wo report in Fig. (3) the lead- 
ing, second-order two-body amplitude for particles of oppo- 
site spins in the degenerate {x,y) orbitals of CUO4. We 
demonstrated [28] that this produces an effective interaction, 
which pushes down the singlet in Eq.(3) and up the triplet by 
|Acu04(4)|. In this way, Acu04(4) can be redefined without 
any reference to the ground state of clusters with a different 
number of particles, and we are free from the objections based 
on a possible Jahn- Teller distortion of odd-iV clusters [33] . 
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FIG. 3. The second order spin-flip diagrams for the two-body 
ampUtude. 



The formation of bound pairs from purely repulsive inter- 
ations was first proposed in a pioneering paper [34] by Kohn 
and Luttinger. They showed that any three-dimensional 
Fermi liquid undergoes a superconducting transition by 
Cooper pairs of very large angular momentum I. A simplified 
view of the Kohn-Luttinger effect is given by considering one 
particle of the pair as an external charge. Then, the screening 
gives rise to a long-range oscillatory potential (Fricdcl oscil- 
lations) due to the singularity of the longitudinal dielectric 
function at 2A:p; here fep is the Fermi wavevector. The strict 
reasoning exploits the fact that the Legendre expansion coef- 
ficients of any regular direct interaction between particles of 
opposite momentum drops off exponentially in I. On the other 
hand, the second-order contribution to the scattering ampli- 
tude falls as l/l'^ and at least for odd I leads to an attractive 
interaction. In the modern renormalization group language 
[35], the second-order correction is obtained by summing up 
the marginal scattering amplitudes of the isotropic Fermi liq- 
uid coming from the so-called Forward channels, including, for 
antiparallel spins, a spin-flip diagram. This scenario does not 
work in the two-dimensional Fermi liquid, but going beyond 
the second-order perturbation theory the Kohn-Luttinger ef- 
fect is recovered [36]. In the Hubbard model, there is first- 
order interaction only for antiparallel spins, but it vanishes for 
W = pairs; in second-order, in the singlet channel, the spin- 
flip diagram is the only one that survives. We found pairing 
in the singlet channel in a variety of models including car- 
bon nanotubes [37] . Hence, the W = pairing mechanism in 



III. MODELLING INTERCELL HOPPING 

We use CUO4 units as nodes of a graph A. The total Hamil- 
tonian is 



^^tot = Ho + H-r, 



(4) 



with 



aGA ia 

i 

where p\ is the creation operator onto the O i = 1,..,4 
of the a-th cell and so on, while Ht is an intercell hopping 
Hamiltonian. The point symmetry group of Ho includes S^^^ , 
with |A| the number of nodes. 

There are many different ways to add an intercell hopping. 
Nevertheless, to preserve the symmetry that produces the 
Acu04(4) < property, we take Hr invariant under the S4 
subgroup of iSl^' . In the following we shall consider a hopping 
term that allows a particle in the i-th O site of the Q-th unit 
to move towards the i-th O site of the /3-th unit with hopping 
integral Tag = jra/jje'^"": 



a,/3eA ia 



(6) 



For A'^ — 2|A| and Tap = 0, the unique ground state consists 
of 2 particles in each CUO4 unit. Section IV (theory) and Sub- 
section VII A (numerical) arc devoted to the intercell hopping 
produced by small {tcpI -C |Acu04(4)|; to study the propaga- 
tion of p pairs we consider a total number of N = 2|A| -|- 2p 
particles. When U /t is such that Acu04 (4) < 0, each pair 
prefers to lie on a single CUO4 and for N = 2|A| + 2p the 
unperturbed ground state is times degenerate (since 

^£ has dimension 2). On the other hand, Subsection VII B re- 
ports the effects of an intercell hopping tcu between Cu sites 
only; this does not break the §4^' symmetry and therefore 
its consequences on pair propagation and SFQ axe drastically 
different. By this sort of models one can study the interaction 
of several pairs in the same system. We are using CUO4 as the 
unit just for the sake of simplicity, but the W = mechanism 
produces bound pairs at different fillings for larger clusters 
[28] too. By replacing CUO4 by larger units one can model 
other ranges of filling fraction. 



IV. LOW-ENERGY EFFECTIVE HAMILTONIAN 
FOR 0-0 INTERCELL HOPPING 

In order to study the propagation of the p added pairs, 
we obtain an effective Hamiltonian by the cell-perturbation 
method with Ho, Eq.(5), the "cell-Hamiltonian" and Ht, 
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Eq.(6), the "intercell perturbation" and by taking into ac- 
count only the low-energy singlet sector. We note that 
the cell-perturbation method was already used in Ref. [38] 
and in Ref. [39] to support the original Anderson's conjec- 
ture [17] on the "low-energy equivalence" between the d — p 
model (proposed by Emery [40]) and the single-band Hub- 
bard model. Despite the analogies with Ref. [39] (like the 
same cell-Hamiltonian and weak O-O links between different 
cells) our intercell perturbation is different and, more impor- 
tant, it is the low-energy sector which differs (one needs to 
consider CUO4 units with 2, 3 and 4 bodies, in contrast with 
0, 1 and 2 bodies of Ref. [39]). 

Let us introduce some useful notation. Let |^g'''(Q)), a = 
1 . . . |A| be the ground state of the q — th CUO4 unit with 
2 < A'' < 4 particles and -E'"-* (A) the corresponding energies. 
Let 5 C A, with |>S| = p, be the set of p CUO4 units occupied 
by four particles and <S its complement in A, that is S = 
A \ <S. When Acu04 (4) < 0, the unperturbed ground state 
with 2|A| + 2p particles can be written as 

and its unperturbed eigenenergy is given by 
Ho\^^) = [pE^°\4) + i\A\-p)E^°\2)] = . 

(8) 

The energies needed to excite the CUO4 cluster from the 
ground state ior U = 5 -{r 6 t [when Acu04(4) « — 0.04t, 
see Fig. (2)1 is w t with 2 bodies and ~ 0.1 t with 4 bodies; 
both are large compared with |Acu04(4)| and no level cross- 
ings take place if \Tafj \ < |Acu04(4)|. 

The perturbation Ht in second order will remove part of 
the 2''x ('^'^ degeneracy. Let |$o) be an exact eigenstate with 
eigenenergy E. Expanding |$o) as 



the fact that the first excited state with 3 particles is w t 
above the ground state for U in the range 5 6 
The energy of the excited states in Eq.(ll) is 



E^°^ = (\S\ - 1)£;^"^(4) + 2£^"^(3) + (|A| - |<S| - 1)£^"^(2) 

and does not depend on the sets T and X>. Within this ap- 
proximation the Schrodinger equation (10) reduces to 



(0), 



7.(0)1 



T.-D S' 



(12) 



where e = E — E^^ and we have disregarded contributions of 
higher order in e. 

The amplitude as = a(ai,...,ap) is totally symmet- 
ric with respect the permutations of the distinct indices 
Qi, . . . ,ap. Letting K.{a) — {(3 € A. : Tap 7^ 0}, after some 
algebra Eq.(12) may be written in the form: 

p 

a(ai, ..,aj-i,l3,aj+i, .., a^)] = £o(ai, . . . , ap). (13) 

This is a Schrodinger equation for p hard-core bosons with a 
complex effective hopping integral; below the J coefficients 
will be calculated analitically and studied as a function of the 
ratio U/t. In Eq.(13), the second term in the l.h.s. describes 
pair propagation, e.g. from unit to an unoccupied unit P; 
in the first term, the system gets back to the initial state after 
virtually exploring unit /3; ni5^j (l ~ '^/3q») takes into account 
that if /3 is one of the occupied units, the particle cannot move 
toward it. 



(9) 



A. Evaluation of the Effective Hopping Integral and 
Selection Rules 



one gets 



{E-Ei°^)as = Y^ 



J2 {i>oU^r\'l\u}i<l>m\IU\<ft ) 



E-E 



(0) 



as> 



(10) 



where {j'l'm)} is a complete set of excited eigenstates of Ho 
and {Em^} their eigenenergies. 

Our crucial approximation is now to truncate the sum over 
the excited states {|3'm.)} considering only the low-energy 

states of the form 



In Eqs.(7)-(11) we omitted the irrep labels of j*,',^^), |*(?^) 
and I in order to avoid a proliferation of indices. Nev- 
ertheless, to calculate the J's we need to reintroduce them. 
The ground state with two particles is non-degenerate; it be- 
longs to the irrep Ai of 54 that coincides with Ai in Cav, and 
the change in notation is: 



-,(2),-4i 



(2),Ai 



(14) 



As far as the ground state with three particles is concerned, 
we recall that it is three times degenerate (apart from the 
trivial spin degeneracy) and belongs to the irrep Ti of 154 that 
in C4v breaks into Bi ® E: 



(r) 



) r = l,2,3 



(15) 



where T C A is the set of |5| — 1 = |T| Cu04 units with 
4 particles, obtained by removing one particle in one of the 
previous |iS| units with 4 particles; in this way we get ]0| = 2 
cells with 3 particles; the remaining |A| — |<S| — 1 s |T U 1*1 
cells have 2 particles. This approximation is legitimated by 



where we may set up our basis such that = Bi , 7^ = 
E^ and T/^' = Ey. 

Finally, the ground state with four particles is two times 
degenerate and belongs to the irrep £ of 154 that in C4v breaks 
into Ai^Bq: 
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(16) 



where f = A\ and E'^'^^ — B2. Useful selection rules may be 
obtained using group theory. By exploiting the invariance of 
Hr under the group 5*4 and omitting, for the sake of clarity, 
the spin indices in the states I'l'o^'' ), it follows that 



(a))i*, 



(3),t; 



('■,3) 



is non- vanishing if and only if 

-'1 ^ -'1 
for = Ai; instead for ft'^ = B2, 

■'1 — ^x, — J^y 



m 

(17) 
(18) 
(19) 
(20) 



In the case s = 1, we have f = A\ and hence by using the 
selection rules in Eq.(18), the matrix element in Eq.(17) con- 
sists of three different contributions coming from the virtual 
3-body states of symmetry B\ ® B\, E^® Ex and Ey ® Ey-. 

r^^[Bi,Bi] = -r.^^(*(3)-^i(a)|pJ„^Jv|>(2''-*Ha)) x 

i 

(*^^J'^M/3)Ip^,.t!*I''-'^(/3)) ; 

Tfp[E„E,] ^ -r„^^{<j'''-(«)|pl,,J*f '^Ha)) x 

i 

with jj, — x,y. Therefore one can solve the Schrodinger equa- 
tion (13) with 



Jc 



a, 13 



: 2 X 



\rf^[BuB,]\^ + \rfp[Ex,Exf + Ir^^l^?., S,]^ 



Acu04(4) 



(21) 



where the factor 2 comes from the spin degeneracy. In Fig. (4) 
we show the trend of J in units of |tP versus U /t. The case 
£;(2) _ ^2 jg gjjjjiiar and may be obtained by group theory. 



J 

r 



U/t 



The above treatment holds for any A; in the next Subsec- 
tion we shall specialize to the case of a one-dimensional chain 
with hopping integrals Tap only between nearest neighbors 
units. 



B. Ring Shaped System 

In order to discuss the propagation of bound pairs and the 
quantization of a magnetic flux, we use a chain with periodic 
boundary conditions and nearest neighbors hopping matrix 
elements Tcr: 



T if /3 = Q + 1, 
r* if /9 = a — 1, 
otherwise 



r = re 



(22) 



For the sake of simplicity we first analize the case with just 
one added pair (p=l) in detail. If the pair belongs to the 
Ai component of the £ irrep the zeroth-order ground-state is 
|A| times degenerate. The effective Schrodinger equation (13) 
reads 



ea{a) = J 2a{a) + e W 't'o a{a + 1) + e 1^1 </>o a(a - 1) 



which is readly solved by Fourier transforming and yields the 
following eigenvalues 



2tt 
W\ 



k + 2- 



k = l, 



,\A\. (23) 



The presence of the factor 2 in front oi (p/cpo implies that the 
model quantizes the fiux in units of /2, like superconducting 
pairs do. Indeed, the ground state energy i?^'''(2|A| -f 2) = 
miufc Sk is strictly periodic in (j) with period (jio/2 for any |Aj > 
2. 

The case |A| = 3 will be used below for a numerical test of 
Eq.(23). The lowest state energies for every quasimomentum 
in the three- unit ring with 8 particles are plotted in Fig. (5) 
versus (j}/4'o- It can be shown, see below, that |A| = 3 is the 
shortest ring showing this effect. 
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FIG. 5. Results of Eq.(13) for the three-unit ring with 
t| = 0.001 t, U = 5 t. Lowest-energy eigenvalues labeled by their 
intercell quasi-momentum are shown versus flux (j>. All energies are 
in t units. 



FIG. 4. J in units of \t\^ versus U/t. The formal divergence as 
(7 ^ is not serious since the effective theory holds for intermediate 
U/t, when |Acu04(4)| is large enough. 
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In the case of p pairs, Eq.(13) yields 



V. SPIN-DISENTANGLED DIAGONALIZATION 



j=i 0=±i i^j 

e V^'^a{ai,..,aj + I3,..,ap)]= sa{ai,. . . ,ap) (24) 

and the model is equivalent to the Heisenberg-Ising spin chain 
governed by the Hamiltonian 



He 



|A| 

E 

a=l 



(25) 

where the a's arc Pauli matrices, spin up represents an empty 
site and spin down represents a pair, q is the so called 
anisotropy parameter and to reproduce Eq.(24) we must 
choose ij = —1. For 77 = 1, we have the isotropic Heisen- 
berg interaction. By performing a Jordan- Wigner transfor- 
mation, the Hamiltonian in Eq.(25) can also be mapped into 
a model of spinloss fcrmions on the ring. In the absence of a 
threading magnetic field {cj> = 0) the problem was originally 
studied by Bloch [41] and then exactly solved by Hulthen [42] 
[in the case = — 1] and Orbach [43] [in the case rj < —1] 
using the Bctlic's hypothesis [44]. A systematic analysis in 
the whole range of parameters was given by Yang and Yang 
in a self-contained series of papers [45]. Here we just recall 
that the model has a gaplcss phase if I?;] < 1, correspond- 
ing to the conducting state, while an insulating phase sets 
in for rj < —1. As in the Id Hubbard model, the "magnetic 
perturbation" (0 ^ 0) does not spoil the integrability and the 
Heisenberg-Ising Hamiltonian remains exactly solvable by the 
Bethe-ansatz method. Let us write an eigenfunction of Hm 
as 



a(ai, Qp) = ^ Ape 



(26) 



where P is a permutation of the integers 1, . . . ,p and Ap are 
p\ coefficients. Shastry and Sutherland [46] have shown that 
the variables kj are given by 



|A|jfc,- = 2nlj + 47r-^ - ^ e{kj,ki) 



with a phase shift 



e{k,q) = 2 tan 



77sin[(fc-g)/2] 



cos[(fc -I- q)/2] - 7?cos[(fc - q)/2] 



(27) 



(28) 



EYom Eqs.(25)-(27) we readly see that the ground state en- 
ergy of the low-energy effective Hamiltonian Hm is periodic 
with period 0o/2, independent of the number of added pairs. 
Moreover Hm is also the appropriate effective model for the 
strong-negative U Hubbard model [47] , for which evidence of 
superconductivity is clear [48]. Thus we conclude that the 
purely repulsive Cu04-Hubbard ring threaded by a magnetic 
field quantizes the flux in a superconducting fashion if the 
number of particles is 2|A| -|- 2p with < p < |A|. 



In order to find the lowest energy eigenvalues and eigen- 

functions when the dimension A/" of the Hilbert space is very 
large one must avoid storing the A/" x A/" Hamiltonian ma- 
trix "Hm-kM- The obvious recipe for sparse matrices, which 
are best handled by the Lanczos method, prescribes storing 
the nonzero elements and their row and column addresses in 
arrays. However for big problems this process is slow, since 
the matrix elements must be referenced to their position in 
'Hj\rxM and must be retrieved each time for later use. One 
can try to improve the situation by projecting on the irreps 
of the symmetry group in order to reduce A/"; however this 
involves building the projected basis and then new Hamilto- 
nian matrix elements, and the new Hamiltonian is much less 
sparse than before, which destroys much of the benefit gained 
by the symmetry. 

Here we exactly diagonalize the |A| = 2 and |A| = 3 ring 
Hamiltonian by the Spin- Disentangled technique, which we 
briefly introduced recently [31], but deserves a fuller illustra- 
tion. This allows to solve the A' x A' many-electron prob- 
lem by storing and handling VJ7 x VJ7 matrices. One of 
its clear advantages is that it does not require the Hamilto- 
nian to be sparse. The method is also much more convenient 
than the usual one, since the x Vf7 matrices are directly 
used by matrix multiplication, which is a fast process, and 
no mapping back onto TiuxM is involved. All the advantages 
of symmetry are gained simply by using a projected starting 
state of the Lanczos chain. To the best of our knowledge, the 
Spin-Disentangled technique was not invented earlier, which 
is somewhat surprising, being a very general method. 

We let A/| + M| = N where is the number of particles 
of spin a; {|<^ao-)} is a real orthonormal basis, that is, each 
vector is a homogeneous polynomial in the p^ and of degree 
Ma acting on the vacuum. We write the ground state wave 
function in the form 



(29) 



a/3 



which shows how the f and [ configurations arc entangled. 
The particles of one spin are treated as the "bath" for those of 
the opposite spin: this form also enters the proof of a famous 

theorem by Lieb [49]. In Eq.(29) Lap is a xmy rectangular 
matrix with m<T = (^'). We let denote the kinetic energy 

vfia X ttIct square matrix of Htot in the basis {\4>aa)}, and nI""^ 
the spin-cr occupation number matrix at site s in the same 
basis {Ns"'' is a symmetric matrix since the \cf>aa)^s are real). 
Then, L is acted upon by the Hamiltonian J?tot according to 
the rule 



U) 



(30) 



In particular for M| = {Sz = sector) it holds K-j = Ki 

(T) fx) 

and Ns = Ns ■ Thus, the action of H is obtained in a spin- 
disentangled way. The generality of the method is not spoiled 
by the fact that it is fastest in the Sz = sector, because it 
is useful provided that the spins are not totally lined up; on 
the other hand, Sz = can always be assumed, as long as the 
Hamiltonian is SU{2) invariant. 



7 



For illustration, consider the Hubbard model with two sites 
a and 6 and two electrons (H2 molecule) each in the 0a or (pb 
orbital. The intersite hopping is t and the on-site repulsion 
U. In the standard method, one sets up basis vectors for the 
Sz = sector 

\tp3) = |06T> <8) |0ai>, \lp4) = |06T> ® |06i), 

One then looks for eigenstates (three singlets and one triplet) 



(31) 



of the Hamiltonian 





t 


t 


°\ 


t 








t 


t 








t 


\ 


t 


t 


uj 



(32) 



Insted of working with 4x4 matrices, we can cope with 2x2 
by the spin-disentangled method using the form in Eq.(29) 
with 



1 




N 



t 

t 



1 



Using Eq.(30), one finds 

Hh,\^}=Y1 Y^iH^^mUlM^lM (33) 



a=o,6 (3=a,h 



with 



Utbi +t(i/'2 + Vs) 



t(-(/)l +^/)4) 
U^i + t{'>l>2 + Vs) 



(34) 



The reader can readily verify that this is the same as applying 
in the form of Eq. (32) to the standard wave function in 
Eq.(31) and then casting the result in the form of Eq.(29). 
Since we can apply /^Hs we can also diagonalize it. 

A/" = 4 for the H2 toy model, but in the 6*2 = sector 
for the |A|=3 ring, M = 1863225 and VJf = 1365, which is 
a clear advantage. Here we have implemented this method 
for the Hubbard Hamiltonian. We emphasize, however, that 
this approach will be generally useful for the many-fermion 
problem, even with a realistic Coulomb interaction, which 
can be suitably discretized. 



A. The practical numerical recipe 

We put a symmetry-adapted trial wave function in the form 
in Eq.(29) and operate the Hamiltonian matrix by Eq.(30); 
each new application introduces a new Lanczos site and we 
can proceed by generating a Lanczos chain. To this end we 
need to orthogonalize to the previous sites by the scalar prod- 
uct given by (^'i|^'2) = Tr(L|L2). In this way we put the 
Hamiltonian matrix in a tri-diagonal form. This method is 



well suited since we are mainly interested in the low-lying 
part of the spectrum. A severe numerical instability sets in 
when the chain exceedes a few tens of sites, i.e. well before 
the Lanczos method converges. Therefore we use repeated 
two-site chains alternated with moderate-size ones. 

In the basis of the sites (of the original cluster) the oc- 
cupation matrices A^i"^' are diagonal with elements equal to 
or 1, simplifying the calculation of the interaction term. 
Moreover, in choosing the trial wave function for small t we 
take full advantage from our knowledge of the ^4 irrep of the 
T = ground state. This speeds the calculation by a factor 
of the order of 2 or 3 compared to a random starting state 
(or even more, if U is large). Typically, starting from a t = 
ground state for the three-unit ring, 24 short Lanczos chains 
were enough to obtain a roughly correct energy and a 20-site 
chain achieved an accurate eigenvalue and an already stabi- 
lized eigenvector. In limiting cases when the results could 
be checked against analytic ones, using double precision rou- 
tines an accuracy better than 12 significant digits was readily 
obtained. On a personal computer with a Celeron CPU the 
three-unit ring with 8 particles required 40 minutes for the 
eigenvalue, and an accurate wave function required less than 
a hour. 



VI. NUMERICAL RESULTS AND DISCUSSION: 
THE TWO-UNIT "RING" 

In this section we analyze the two-unit ring with 

Acu04 (4) < 0. Here a complex r is equivalent to a real one, 
I. e. no magnetic field can thread the system. This is peculiar 
of the ring with two clusters, since each CUO4 is on the right 
as well as on the left of the other. Indeed, we have 

Hr = iT + t") ^ {plj^PiJa +P2,jaPl.J>) = 

= 2|t| cos ( ^ j ^ {p\j^P2,ja + plj^PlJa) (35) 

hence, the ground state energy J5'-°-'(0, |r|) as a function of 
the flux d) and of the modulus of r satisfies 



E^°\4>,\t\) = E^">{0, |t|cos^ 



(0)/ 



(36) 



E^'^\(f),\T\) has a local maximum at = 0o/2 (a property 
which is independent of the number of particles in the two- 
unit ring) since Hr in Eq.(35) vanishes; see Fig. (6) for the 
case of 2|A| -(-2 = 6 particles. 

The translational invariance allows to label any state by 
the crystal momentum p = nhk, k = 0,1. We introduce the 
short-hand notation |-|-) and |— ) for the components of the 6- 
body non-interacting ground-state multiplet with A; = and 
k = 1 respectively 



l±> = 



1 

71 



*(*)(l))|*^^^(2))±|*^^^(l))|*r(2)> 



-,(2) 



(2)/ 



-,(4)/ 



We have 



(*(^)(i)K*r(2)|J^.|->=o 



>(3) 



(37) 



(38) 
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(the singlet projection of |*^^'(1))|*^^^(2)) has fc = quan- 
tum number) and there is no second order correction in the 
k — 1 subspace as shown in Fig. (6). In the fc = subspace 
the correction is proportional to |rp/Acu04 for small \t\, 
in agreement with the analytical predictions. As shown in 
Fig. (6), the maximum at = (j}o/2 is not a cusp, as there no 
level crossing is found. 



In Fig. (7) we report the numerically exact results for a ring 
with three clusters, with U — 5t and |r| in the range from 
to 0.006t. The ground state energy and the first excited level 
depend quadratically on |t| for small |r|. As expected from 
Eq.(23), the lower eigenspace has fe = 0, while the first excited 
one contains the states with fe = 1 and 2. Differently from the 
two-unit ring, the first excited level receives a second-order 
correction. 
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FIG. 6. Energy of the ground state (fe = 0) and of the first 
excited state (fc = 1) with 2|A| + 2 = 6 particles as a function of 
0/</>O. Here U = 5t, (Acu04(4) ~ -0.04258 1) and |r| = 0.001 i. 
The energies are in units of t. 



-8.9882 



-8.9883 











































>< 


























k=0 








/'\ 


k=l 








=2 

















0.5 1 

FIG. 8. Numerical results for the low energy states of the 
three-unit ring , as a function of the concatenated magnetic flux. 
Here U = 5 t, (Acu04(4) ~ -0.04258 i) |t| = 0.001 1. The energy 
is in units of t. 



VII. NUMERICAL RESULTS AND DISCUSSION: 
THE THREE-UNIT RING 

A. O-O intercell hopping 

In this section we consider the three-unit ring focusing the 
attention on the case Acu04 (4) < and total number of par- 
ticles 2| A| +2 = 8. The switching on of the hopping r between 
the O sites breaks the symmetry group Csv <8> into Csv <8> S4 
for real r; in a magnetic field (complex r), this further breaks 
into C3 ® 54. Real r lifts the degeneracy between the fc = 
subspace and the subspaces fe = 1 and fe = 2 of C3 (as usual fe 
is related to the crystal momentum p = 2nhk/3 in this case), 
but cannot split fe = 1 and 2 because they belong to the 
degenerate irrep of Ca^; complex t resolves this degeneracy. 
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FIG. 7. Energies of the ground state (fc = 0) and of the first 
excited state (fc' = 1 and 2) of the three-unit ring with 8 particles 
versus |t| for U = bt. The energies are in units of t, the parameter 
|t| is in units of W~^t. 



The three-unit ring is the smallest ring where we can insert 
a magnetic flux by r = l^le*^, 9 = ^{(l)/(j)o)- The energies 
of the three ground-state multiplet components are reported 
in Fig.(8) for |r| < |Acu04(4)| and U = 5t. At (f) = the 
ground state belongs to the fc = subspace, while the first 
excited levels have fe = 1 and 2. Their spatial degeneration 
is fully lifted: the fc = 1 level increases while the k = 2 level 
decreases up to = 0o/2. As (j) increases, the ground state en- 
ergy grows quadratically in cj> (diamagnetic behaviour). Near 
= 0o/4 we find a level crossing between fe = and fe = 2, 
while &t (j> = (j>o/2, k = Q becomes degenerate with fc = 1 and 
the ground state energy is in a new minimum belonging to the 
k = 2 subspace: a sort of "restoring" of the <j) = Q situation is 
taking place as in the BCS theory [50]. Indeed, at <;!> = 4>o/2 
the symmetry group is C31, ® S4 where Cs^ is isomorphous 
to Csv (reflections a are replaced by ag, where g is a suit- 
able gauge transformation). This feature was also found in 
other geometries [18] [32]. In the region from (f) = <po/2 and 
= 00 we numerically verified that Ek-2{<i>) = -Efe=2(0o — 0), 
Ek=o{(j)) = Ek=i{<j)o-<f}) and Ek=i{(f>) = Ek=o{4>o - 4>) , where 
Ek{(l>) is the ground state energy in the fe sector. The compari- 
son of the numerical results shown in Fig. (8) with the analytic 
ones in Fig. (5) supports the accuracy of the cell-perturbation 
scheme proposed in Section IV. Thus, the dressed W = 
pair screens the vector potential as a particle with an effec- 
tive charge e* = 2e does. At both minima of E'''^\(f)) we have 
computed A3_unit(8) ~ — 10^'^t. Here, the half-integer AB 
effect is actually SFQ. 

Fulfilling the conditions Acu04(4) < and |r| ^ 
|Acu04(4)|, we varied U and jrj and found analogous trends 
for the ground state energy. Increasing jrj with fixed 
Acu04 (4) lowers the central minumum and depresses the two 
maxima. On the other hand, if |Acu04(4)| decreases at fixed 
|t| the central minimum and the side peaks are affected in a 
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similar way. This is reasonable since the perturbative param- 
eter is |t|/|Acu04(4)|. 

The three-unit ring also enables us to study persistent dia- 
magnetic currents carried by bound pairs screening the mag- 
netic flux. We calculated the expectation value for each k 
of the total current operator as a function of the flux. The 
current operator [51] 



The absence of SFQ in Fig. (10) is a further evidence of the 
repulsion-driven pairing mechanism discussed in Section II. 



dHu 



d(j} h\A\ 



* t 

Pa.ic 



Pa + l.icT ) 



(39) 



yields a gauge invariant average I. By expanding I in Eq.(39) 
in powers of (j) near (j> ~ one may identify the paramagnetic 
and the diamagnetic contributions with the zeroth and the 
first order terms respectively [52] . The results are reported in 
Fig. (9); the current is proportional to the flux derivative of the 
ground-state energy [see Fig. (8)] according to the Hellmann- 
Feynman theorem. Near (j) = the system generates a dia- 
magnetic current which screens the threaded magnetic field. 
When <f) exceedes a critical value ~ (po /4, a breakdown of the 
ground state occurs. This corresponds to a discontinuity of 
the current which changes sign; then the current enhances the 
external field. At — (f)o/2 the current vanishes again. In- 
deed, like at = 0, the eigenfuctions may be choosen real [H-r 
at = (/)o/2 is obtained from H-r at = by reversing the sign 
of four 0-0 bonds connecting two nearest neighbours units]. 
Thus, near (j>o /2 the magnetic flux is still a small perturbation 
with respect to a new real intercell hopping Hamiltonian and 
the current correctly screens the new magnetic field. From 
Fig. (9) we see that the maximum value of the diamagnetic 
current is of the order of 1 10 nano Ampere if t = 1 eV and 
the ratio I /{(j)/cj>o) ~ e|r|//i near = 0. 
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FIG. 9. Total current for the three-Cu04 ring, as a function of 
the magnetic flux. Here U = 5 t, \t\ = 0.001 1. The current is in 
units of e\T\/h. The thick fine marks the ground state current. 



In Fig. (10) we show the trend of the ground state energy 
in each k sector for the non-interacting {U — 0) three-unit 
ring. In this case there is no pairing in Cu04 and indeed 
the ground state energy is linear in the field at small fields 
(normal Zeeman effect). The lowest state is A; = 2 through- 
out. Interestingly, the three-unit ring concatenated with half a 
flux quantum would be diamagnetic, but the absence of a sec- 
ond minimum shows that it would be Larmor diamagnetism. 
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FIG. 10. Low energy states of the three-unit ring versus 
Here (7 = and |t| = 0.001 1. The energy is in units of t. 



B. Cu-Cu intercell hopping 

We can alternatively model the three-unit ring by connect- 
ing only the central Cu sites of the constituent Cu04 units 
with a hopping term rcu; in order to study the propagation 
of a bound pair we again assume the total number of parti- 
cles 2|A| 4-2 = 8. The full system threaded by the flux has 
a. Cs ® S'i symmetry because the O sites are not involved 
in the intercell Hamiltonian. Again, we consider the case 
Acu04(4) < and ask if now the W = G bound pair can 
screen out the magnetic flux. 
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FIG. 11. Ground state energy E of the three-unit ring in 
units of t, as a function of the concatenated magnetic flux. E is 
^-independent (see text). Here U = 5 t, I'^Cu 

= O.lt. 



We find that rcu produces much smaller effects than r; 
for |tcu| <C |Acu04(4)| the energy eigenvalues are outright 
(/> independent to great accuracy. Therefore we considered 
TCul = O.lt; still the dependence of the ground state energy 
is weak, see Fig. (11). At = 0, the correction [53] due to [rcuj 
to ground state energy is ~ 10~^t. Moreover, remarkably, the 
system behaves as a paramagnet. 

The reason of this unusual behavior is the following. There 
is no flux-induced splitting of the three k levels because the 
W = pair is strictly localized by the local symmetry. Indeed 
the 54 label of each Cu04 unit is a good quantum number. 
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No SFQ is observed because the screening of the magnetic 
field by the bound pair is forbidden. The small correction 
to the ground state energy comes from a second-order pro- 
cess. Starting e.g. with an unperturbed state |4, 2, 2), in 
which the bound pair is localized on the first cluster, the cor- 
rection involves virtual states |3, 3, 2) and |3, 2, 3) in which a 
totalsymmetric particle jumps forth and back on the nearby 
clusters. This process occurs with a small amplitude because 
of a severe energy misfit. This is particularly clear at weak 
coupling, when the lowest-energy Ai particle of the first clus- 
ter must hop to antibonding Ai orbitals of the nearby clus- 
ters; the amplitude of this process is further reduced by the 
overlap of these orbitals with the localized Cu one. How- 
ever, such virtual processes are insensitive to the flux. Any 4> 
dependence arises from third-order corrections (order |A| in 
general). Indeed, the A-i particle must go virtually around 
the trip clockwise or anticlockwise. In the ground state, of 
course, it chooses the wise in such a way to gain energy from 
the magnetic field. This is why a paramagnetic dependence 
on the flux is seen in Fig. (11) and the correction goes like 
—(jP' at small 4>. This is interesting because it shows how the 
local symmetry can hinder the tunneling of bound pairs car- 
rying conserved quantum numbers; SFQ is not a necessary 
consequence of superconductivity if the pairs are not total- 
symmetric. 



VIII. CONCLUSIONS 

We propose a Hubbard model with on-site repulsion de- 
flned on a graph A with 5-site C4«-symmetric clusters as nodes 
in order to study the response of the system to a threading 
magnetic flux. For ring-shaped systems and weak O-O links 
we find a half-integer AB effect which is unambiguosly inter- 
preted as SFQ. The key ingredient is the W = Q pair which is a 
two-body singlet cigenstate of HcuOi without double occupa- 
tion and is formed by mixing degenerate one-body states. In 
the interacting problem the W = pair becomes a bound pair 
when four particles lie in the Cu04 cluster. The pairing mech- 
anism is due to an effective attractive interaction mediated by 
repeated electron-hole exchanges with the Fermi sea. Thus, 
SFQ may be found in purely repulsive Id Hubbard models if 
the nodes are represented by a non-trivial basis. Focussing on 
the low-energy sector, we find a simplified description of the 
model in terms of an effective hard-core boson Hamiltonian 
that can be solved exactly for ring-shaped systems and arbi- 
trary filling. Further, we show that the boson Hamiltonian is 
equivalent to the well-known [42] Heisenberg-Ising spin chain 
with an antiferromagnetic anisotropy parameter 77 = —1. 

The analytic results are well confirmed by the numerical 
findings for the two- and three-unit ring (14,400 and 1,863,225 
configurations). To this end, we have recently introduced a 
new exact-diagonalization technique. By disentangling spin- 
up and spin-down we reduce the size of the matrices that 
must be handled considerably; for spin-unpolarized systems 
the matrix dimension is the square root of the overall size of 
the Hilbert space. The method allows large further reduc- 
tions by exploiting the symmetry; this can be done in several 
ways. One can first diagonalize Dirac's characters and then 
apply the spin-disentangled technique with a smaller function 



space; alternatively, one can set up the spin-disentangled tech- 
nique by a variational approach, a projection operator (i/") 
or Lanczos method, starting with a trial state belonging to a 
well defined symmetry. 

We have computed the ground state energy and the induced 
supercurrcnt as a function of the trapped flux in the case of 
weak 0-0 links. We have also studied the effect of direct 
intercell Cu-Cu links; in this case, bound pair propagation 
is hindered by symmetry, because each unit must keep its 
own S4 irrep. Hence, the unusual situation arises when the 
threading flux is not screened by the superconducting pairs 
and a paramagnetic response prevails. 
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